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Abstract.  There  is  strong  physical  evidence  that  a  full  treatment  of  differential-algebraic 
equations  should  incorporate  solutions  with  jump  discontinuities.  It  is  shown  here  that  for 
semilinear  problems  the  setting  of  distributions  allows  for  the  development  of  a  theory  where 
indeed  such  discontinuities  may  occur.  This  approach  also  settles  the  problem  of  inconsistent 
initial  conditions  in  a  very  simple  way.  On  the  other  hand,  new  issues  arise  as  not  only 
uniqueness,  but  even  countability  of  the  number  of  solutions  of  Initial  value  problems  may 
now  be  lost.  A  physically  motivated  but  purely  mathematical  selection  procedure  to  overcome 
this  difficulty  is  discussed  in  Part  11  of  this  paper. 


1.  Introduction. 

The  investigation  of  quasilinear  differential-algebraic  equations  (DAE’s)  in  R", 

(1.1) 

given  in  [RRh2]  and  [RRh3]  (see  also  [CD],  [SDe]  and  [T])  has  revealed  the  presence  of 
singularities,  notably  impasse  points,  be3^ond  which  classical  solutions  cannot  be  continued. 
Thus,  under  the  assumption  that  the  system  (1.1)  governs  the  evolution  of  the  state  variable 
X  at  all  times,  discontinuous  and  hence  nonclassical  solutions  of  (1.1)  may  play  a  key  role 
in  such  problems. 

For  linear  time-dependent  problems 

(1.2)  A{i)x  +  B{t)x  =  h 
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where  b  G  {T>'(J))^,J  C  R  is  an  open  interval,  and  x  denotes  the  derivative  of  x  in  the 
sense  of  distributions,  a  rather  complete  theory  can  be  developed  ([RRh4]).  More  specifi¬ 
cally.  when  b  is  “almost”  a  function,  initial  value  problems  associated  with  (1-2)  continue  to 
make  sense  and  have  a  unique  solution  ([RRh5]).  The  admissible  class  of  distributions  for 
the  validitv  of  such  results  includes  functions  with  jump  discontinuities.  It  turns  out  that 
for  index-one  problems,  the  solutions  are  also  functions  with  jump  discontinuities  (i.e.,  not 
classical  ones)  while  for  higher  index  problems  they  may  also  exhibit  an  “impulsive”  part; 
that  is,  contain  a  linear  combination  of  derivatives  of  Dirac  delta  distributions.  It  is  note¬ 
worthy  that  the  jumps  in  the  solutions  me  calculable,  and  hence  that  these  discontinuous 
solutions  are  characterized  as  completely  as  the  classical  ones. 

It  becomes  evident  that  an  attempt  at  resolving  the  ambiguity  created  by  the  presence 
of  impasse  points  in  (1-1)  should  be  made  via  the  concept  of  distribution  solution.  This  is 
the  topic  of  this  paper  for  the  special  case  of  semilinear  DAE’s  of  index  one 

(1.3)  Ait)i  ^G{t,x). 

The  restriction  to  this  case  is  justified  by  the  following  consideration;  If  -4  is  sufficiently 
smooth,  the  product  A(t)x  is  well  defined  for  any  distribution  x,  whereas  an  expression 
such  as  A(x)i  usually  makes  no  sense  even  when  x  is  a  piecewise  smooth  function  with  a 
jump  discontinuity.  Indeed,  the  discontinuity  in  x  induces  a  discontinuity  in  A(.t)  which 
therefore  cannot  be  multiplied  by  the  Dirac  delta  distribution  in  x  arising  from  the  jump 
of  X.  There  are  ways  to  circumvent  (in  part)  this  difficulty,  but  only  at  the  expense  of 
additional  technicalities  that  we  prefer  not  to  consider  at  this  time. 

Fortunately,  the  restriction  to  semilinear  problems  does  not  significantly  affect  the  range 
of  applications,  as  many  concrete  problems  with  impasse  points  or  other  singularities  even 
involve  a  constant  matrix  A  in  (1.3).  But  it  is  important  that  the  problem  be  of  index 
one,  for  the  semilinear  structure  is  not  preserved  by  the  reduction  procedures  involved  in 
higher  index  problems. 

Remark  1.1:  In  fact,  no  reduction  procedure  described  in  the  literature  applies  without 
restriction  to  distribution  solutions,  except  that  given  in  [RRh4]  for  the  linear  case.  The 
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other  general  techniques  (see  [RRhl],  or  [CG]  developing  another  point  of  view  based  on 
augmented  systems)  implicitly  use  the  fact  that  the  solutions  and  their  derivatives  have 
pointwise  values,  and  hence  are  not  applicable  to  distributions  as  they  stand.  □ 

Usually,  nonautonomous  problems  such  as  (1.3)  are  more  conveniently  handled  within 
the  framework  of  autonomous  systems  by  changing  the  variable  x  into  (t,  x)  and  adding  the 
equation  t  =  1.  But  because  this  transforms  a  semilineax  problem  into  a  quasilinear  one, 
it  is  no  longer  appropriate  to  do  so  when  solutions  are  sought  in  the  sense  of  distributions. 
Since  this  work  relies  rather  heavily  upon  prior  results  established  for  classical  solutions 
and  hence  developed  mostly  for  autonomous  problems  (for  the  reason  just  mentioned),  it 
is  important  to  rephrase  these  results  for  the  case  of  the  nonautonomous  system  (1.3),  thus 
making  it  explicit  how  the  variable  t  is  involved  in  various  places,  notably  in  the  definition 
of  impasse  points.  This  is  done  in  the  next  section,  where  some  additional  properties 
are  eJso  highlighted,  as  they  become  important  later  (but  were  not  involved,  hence  not 
emphasized,  in  the  treatment  of  classical  solutions). 


Section  3  deals  with  the  discontinuous  solutions  of  (1.3).  As  expected,  the  setting 
of  distributions  allows  the  solutions  to  jump  at  impasse  (and  other)  points.  But  this 
pleasant  conclusion  comes  along  with  another,  much  less  welcome  one:  Once  the  setting  of 
distributions  has  been  adopted  and  discontinuities  have  become  possible,  they  abuse  the 
opportunity.  Jumps  may  eilso  occur  in  places  where  there  was  no  reason  to  expect  them 
before,  and  simple  examples  show  that,  now,  initial  value  problems  may  have  uncountably 
many  solutions.  It  is  the  purpose  of  Part  II  of  this  paper  to  show  how  this  difficulty  may 
be  overcome.  Meanwhile,  in  Section  4  it  is  shown  that  distribution  solutions  also  allow  for 
a  simple  answer  to  the  well  known  question  regarding  “inconsistent”  initial  values. 


Remark  1.2:  Any  discussion  of  the  discontinuous  solutions  of  (1.3)  must  include  jumps 
at  singularities  since  this  is  where  jumps  are  widely  known  to  occur  usually.  Then  the 
justification  of  the  distribution  approach  requires  precise  information  about  the  behavior 
of  the  classical  solutions  of  (1.3)  in  the  vicinity  of  singularities.  For  most  singularities, 
this  information  is  still  lacking,  but  in  the  case  of  impasse  points  it  is  hidden  in  the 
(combined)  works  [R]  and  [RRh2].  This  is  why  our  discussion  of  jumps  at  singularities  is 
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limited  to  impasse  points,  and  presumably  one  of  the  reasons  why  the  very  natural  idea 
of  introducing  discontinuous  solutions  of  (1.3)  via  distributions  has  not  appeared  much 
earlier  in  the  literature.  □ 

The  “glut”  of  distribution  solutions  appears  to  be  a  serious  shortcoming  of  the  distribu¬ 
tion  approach,  which  has  no  analog  for  linear  problems.  In  the  second  part  of  this  paper, 
we  try  to  resolve  this  shortcoming,  without,  of  course,  repudiating  the  concept.  Evidently, 
depending  upon  the  “physical”  setting,  many  of  the  unwanted  distribution  solutions  may 
be  ruled  out  for  lack  of  relevance,  but  such  a  decision  is  based  upon  a  non-mathematical 
argument.  In  Part  11  of  the  paper,  we  will  show  that  a  selection  can  be  made  accord¬ 
ing  to  physicedly  motivated  but  purely  mathematical  criteria.  The  idea  is  simply  that 
acceptable  solutions  must  be  consistent  (in  a  way  to  be  defined)  with  a  given  class  of  per¬ 
turbations.  In  practice,  the  admissible  perturbations  are  dictated  by  the  physical  origin 
of  the  problem.  This  concept  of  consistency,  called  “"P-consistency”,  where  V  stands  for 
“perturbation”,  makes  standard  and  novel  connections  with  singular  perturbation  theory. 
Standard  because  the  underlying  idea  is  that  solutions  of  the  unperturbed  problem  should 
be  “approximated”,  locally  at  least,  by  solutions  of  the  perturbed  problems,  and  novel 
because  the  known  criteria  for  this  property  to  be  true  are  used  to  sort  out  the  solutions 
of  the  imperturbed  problem,  and  to  discard  many  of  the  spurious  distribution  solutions. 

Naturally,  almost  all  work  devoted  to  the  discontinuous  solutions  of  (1-3)  stresses  con¬ 
nections  with  singular  perturbation  theory  (with  various  degrees  of  emphasis).  However, 
to  summarize,  the  general  trend  has  been  to  define  the  discontinuous  solutions  of  (1.3)  as 
pointwise  or  other  limits  of  solutions  of  ODE  perturbations  (see  e.g.  [SDe]).  In  sharp  con¬ 
trast,  the  distribution  approach  permits  to  consider  such  discontinuous  solutions  without 
any  reference  to  perturbations.  Perturbations  become  involved  when  it  comes  to  select¬ 
ing  the  meaningful  solutions,  but  then  the  approximation  criterion,  difficult  to  check  in 
practice,  may  be  replaced  by  a  weaker  and  much  more  convenient  eigenvalue  condition. 
Furthermore,  DAE  rather  than  ODE  perturbations  can  be  used  for  the  selection  procedure. 
This  is  especially  useful,  as  a  number  of  physically  motivated  perturbations  arise  as  DAE’s 
cind  not  ode’s. 
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Expanded  examples  (from  electrical  network  theory)  are  given  at  the  end  of  Part  II. 
Since  these  examples  also  illustrate  various  points  made  in  this  first  part,  we  have  not 
included  here  any  other  (physically  motivated)  examples. 


2.  Geometrically  Nonsingular  Serailinear  DAE’s  of  Index  One. 

The  material  presented  in  this  section  is  in  part  condensed  from  the  articles  [RRhl]  and 
[RRh2]  and  specialized  to  the  case  of  semilinear  DAE’s;  that  is,  to  problems  of  the  form 

(2.1)  A{t)x  —  G{t,x), 

where  A  €  C°°(  J;£(R”)),  G  €  C^{J  xR";R")  and  J  C  R  is  an  open  interval.  Consistent 
with  the  requirements  that  (2.1)  is  a  DAE  and  not  an  explicit  ODE,  we  shall  assume 
throughout  that 

(2.2)  rank  A(<)  =  r  ,  0  <  r  <  n,  V#  €  /• 

Furthermore,  since  all  the  results  in  this  paper  involve  only  loced  assumptions  in  the  “time” 
variable  f,  it  is  not  restrictive  to  suppose  that  the  interval  J  has  been  shrunk  so  that  there 
exists  a  common  complement  Z  of  dimension  n  —  r  to  all  the  spaces  rge  A(t),t  G  J: 

(2.3)  R”  =  rge  A{t)  ©  Z,  Vf  €  J. 

We  shall  call  Q{t)  6  £(R”,Z)  the  projection  onto  Z  associated  with  the  decomposition 

(2.3) .  By  elementary  arguments  it  follows  that 

(2.4)  QeC'“(J;£(R”,Z)). 


As  is  well-known,  condition  (2.2)  alone  is  not  sufficient  to  provide  a  satisfactory  existence 
theory  for  the  classical  solutions  of  the  D.4E  (2.1).  Accordingly,  we  assume  also  that 

(  The  mapping 


(2.5) 


(f,x)  e  J  X  R"  : — >  (5(f)G(<,x)  G  z 

is  a  submersion  at  each  point  of  its  zero  set. 
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Condition  (2.5)  is  readily  seen  to  be  independent  of  the  space  Z  in  (2.3),  and  it  implies  at 
once  that  the  set 


(2.6)  IV  =  {(t,  a:)  G  J  X  R"  :  Q{t)G{t,  a:)  =  0}  =  {(t,  x)  €  J  x  R"  :  G{t,  x)  G  rge  A{t)] 

is  a  closed  (r  +  l)-dimensional  C°°  submanifold  of  /  x  R"  (also  independent  of  Z).  Now, 
the  mapping 

(2.7)  {t,  x,p)  G  J  X  R"  X  R"  I — >  A(t)p  ~  G{t,  x)  G  R”, 


with  the  derivative 

(2.8)  (r,  /i,  g)  G  R  X  R”  X  R"  i — >  rDtA{t)p  +  A(t)g  —  TDtG{t,  x)  —  DxG(t,  x)h  G  R”. 

is  a  submersion  at  each  point  {t,x,p)  of  its  zero  set.  Indeed,  we  have  {t,x)  G  W  and 
hence  for  every  u  G  R”  there  is  (by  (2.5))  a  (r, /j)  G  R  x  R"  such  that  TDtQ{t)G(t,x)  + 
TQ(t)DtG{t,x)  +  Q(t)DxG{t,x)h  =  ~Q(t)u.  From  the  identity  (5(t)a4(t)  =  0,  we  infer 
that  DtQ{t)A{t)  —  -Q{t)DtA{t).  Together  with  the  relation  G(t,  x)  =  A(t)p,  this  yields 
DtQ{t)G{f,x)  =  —Q{t)DtA{i)p,  and  hence 

u  -  TDtA{t)p  +  rDtG(t,  x)  +  DxG{t,  x)h  G  rge  A{t), 

i.e.  there  is  g  €  R"  such  that 

TDtA{t)p  +  A(t)q  -  rDtG{t,  x)  —  DxG{t,x)h  =  u. 

Since  u  G  R"  is  arbitrary,  it  follows  that  (2.8)  is  surjective.  (The  converse  is  also  true: 
If  (2.8)  is  surjective  at  the  points  of  the  zero  set  of  (2.7),  then  (2.5)  holds.  The  proof  is 
trivial  modulo  the  remark  that  DiQ{t)G(t,x)  —  —Q(t)DtA(t)p  already  used  before.) 

It  follows  from  all  this  that  the  set 

{(t,  r,  p)  G  J  X  R”  x  R”  :  A{t)p  —  G{t,  i)  =  0) 
is  a  closed  (n  +  l)'dimensional  C°°-submanifold  of  J  x  R"  x  R",  and  hence  that  the  set 

(2.9)  M  =  {(f,  X,  l,p)  G  J  X  R”  x  R  X  R"  :  A{t)p  ~  G{t,  x)  =  0} 
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is  a  closed  (n  +  l)'dimeiisional  C°°  submanifold  of  J  x  M"  x  R  x  R".  Identifying  J  x  R”  x 
R  X  R"  ~  r(  J  X  R”)  (tangent  bundle  of  the  first  two  factors)  we  see  that  M  can  be  viewed 
as  a  submanifold  of  T(J  x  R").  On  the  other  hand,  from  the  embedding  W  C  J  x  R" 
we  infer  that  TW  C  T(  J  x  R"),  and  the  set-theoretic  intersection  TW  0  M  makes  sense. 
Note  also  that 

(2.10)  W  =  n(M), 
where  11 :  r(J  x  R”)  — >  J  x  R”  is  the  canonical  projection. 

Definition  2.1.  The  pair  {t,x)  €  J  x  R”  is  consistent  with  the  DAE  (2.1)  if  (t,x)  G 
Tl{TW  n  M)  C  W  (see  (2.10)).  Accordingly, 

(2.11) 

is  called  the  set  of  consistent  points  for  the  DAE  (2.1). 

The  third  and  final  assumption  needed  to  obtain  an  existence  and  uniqueness  theory 
for  the  DAE  (2.1)  is 

(2.12)  (t,r)  €  Q{t)D,G{t,x)\u..Mo  ^  GI(ker  A(t),Z). 

It  is  once  again  straightforwcird  to  check  that  condition  (2.12)  is  independent  of  the  choice 
of  the  space  Z  in  (2.3).  A  useful,  equivalent  formulation  of  condition  (2.12)  is  contained 
in  the  following  proposition. 

Proposition  2.1.  We  have 

(2.13)  ir  :=  {(t,x)  G  W  :  g(t)D,(?(i,  €  GLiker  A{t),Z)}  C  W^ 

In  particular,  condition  (2.12)  holds  if  and  only  if 

(2.14)  W=  =  W". 

As  a  result,  if  (2.12)  holds,  then  is  an  open  subset  ofW. 

Proof.  Let  (<,x)  €  W*^  be  given.  Since  C  W,  there  exists  a  p  G  R”  such  that 

A(t)p  =  (j(t,x).  Because  QG  takes  values  in  Z,  we  have  Dt(QG)(t,x)  €  Z,  and  by 
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the  surjectivity  of  Q{t)DxG(t,  there  is  a  fc  €  her  A{t)  such  that  Q{t)DxG{t,  x)k  — 

-Q(t)DxG{t,x)p-  Dt{QG){t,xy,  that  is, 

{2.15)  Dt{QG){t,  x)  +  Q{t)DxG{t,  x){p  +  k)  =  0. 

Now,  from  (2.5)  and  (2,6), 

T(t,x)W  =  keTD(QG){t,x)  =  {(r,h)  e  M  X  M"  :  TDt{QG){t,x)  +  Q{t)DxGit,x)h  =  0}, 

so  that  relation  (2,15)  also  reads  (l,p  +  k)  £  T(^t,x)^  hence  (t,a;,  l,p  +  k)  £  TW . 
Also,  (t,  X,  1,  p)  e  M  since  A{t){p  +  k)-G{t,x)  =  A{t)p  -  G{t,  x)  =  0  (see  (2.9)  and  recall 
k  £  ker  A(t)).  Thus,  {t,x,  l,p  +  k)  £  TW  n  M  and  therefore  (t,x)  £  Il{TW  H  M)  =  VT'". 

Since  (2.12)  amounts  to  the  inclusion  C  the  above  proves  the  equivalence  of 
(2.12)  and  (2.14).  An  elementary  contradiction  argument  shows  that  is  open  inW .  □ 

Definition  2.2.  The  semilinear  DAE  (2.1 )  is  geometrically  nonsingular  of  index  1  if  the 
conditions  (2.2),  (2.5)  and  (2.12)  (or,  equivalently,  (2.14))  hold. 

The  relex^nce  of  these  concepts  to  the  existence  of  classical  solutions  of  the  DAE  (2.1) 
is  provided  by  Theorem  2.1  below.  The  given  proof  only  establishes  the  connection  with 
more  general  results  in  [RRhlj  or  [RRh2]  from  which  it  can  be  derived. 

Theorem  2.1.  Let  the  DAE  (2.1 )  be  geometrically  nonsingular  of  index  1.  Then: 

(i)  If  X  £  (/;  R")  is  a  solution  of  (2. 1 ),  we  have  (t,  x(t))  £  W‘^,  for  alii  £  J. 

(ii)  Conversely,  given  (to,  2:0)  €  W  and  after  shrinking  J  about  to  if  necessary,  there  is  a 
unique  solution  x  6  C^(  J;R"),  actually  of  class  C°°,  such  that  x(to)  —  •'I’o- 

Proof.  By  adding  the  equation  t  =  1  and  setting  x  =  (t,x),  the  D.4.E  (2.1)  is  transformed 
into  the  quasilinear  DAE 

(2.16)  i(x)x  =  G(x), 

where 

lit))  ’ 
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and  the  initiaJ  condition  x{to)  =  xq  becomes 


(2.18)  x(to)  -  {to,xo). 

Conditions  (2.2),  (2.5)  and  (2.12)  amount  to  saying  that  the  autonomous  quasilinear 
DAE  (2.16)  is  geometrically  nonsingular  of  index  1  in  the  sense  of  [RRh2].  Here  it  is  worth 
mentioning  that  condition  (2.12)  is  equivalent  with 

z  €  =>  rank  =  r  +  1(=  rank  A(z)), 

which  follows  from  the  characterization  TiW  =  kerZ)(QG)(z)  and  the  remark  that  (2.12) 
amounts  to 

X  kerA(z)|j,.^^  =  {0}. 

Now  the  result  is  a  direct  consequence  of  [RRh2,  Theorem  4.1],  except  for  the  C°°  - 
smoothness  of  the  solution  which  is  obvious  from  the  proof  given  there  when  A  and  G  are 
of  class  C°°.  Alternatively,  the  theorem  also  follows  from  [RRhl,  Theorem  6.1]  (in  [RRhl], 
geometrically  nonsingular  DAE’s  are  simply  called  “nonsingular”).  □ 

The  terminology  “geometrically  nonsingular”  in  Definition  2.2  is  justified  by  the  fact 
that  the  various  sets  carrying  the  information  needed  in  the  existence  theory  are  equipped 
with  a  “natural”  differentiable  structure  (the  proof  of  Theorem  2.1  eventually  relies  upon 
standard  theory  of  ODE’s  on  manifolds).  Thus,  the  problem  does  not  exhibit  any  visi¬ 
ble  (geometric)  singularity.  But  Definition  2.2  also  allows  for  invisible  (or  algebraic;  see 
[RRh2])  singularities,  as  we  now  explain. 

Suppose  that  the  DAE  (2.1)  is  geometrically  nonsingular  of  index  1.  From  Theorem 
2.1,  no  path  (<,z(<))  goes  through  a  point  of  the  set  W  \  (closed  in  W),  if  z  is  a 
solution  of  (2.1).  But  such  points  may  well  lie  at  the  “beginning”  or  the  “end”  of  such 
a  path  (t,  x{t)),  for  unlike  in  explicit  ODE  theory,  here  trajectories  may  stop  abruptly  at 
points  reached  in  finite  time.  Impasse  points^  as  defined  below  (Definition  2.3)  are  the 
most  frequently  encountered  points  of  this  type. 

Impasse  points  are  defined  in  [RRh2]  for  general  quasilinear  DAE’s.  Hence,  in  principle, 
a  definition  for  an  impasse  point  of  (2.1)  can  be  obtained  by  applying  that  definition  to  the 
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DAE  (2.16)  equivalent  to  (2.1).  However,  the  general  definition  of  impasse  points  makes 
reference  to  the  rather  unintuitive  concept  of  the  intrinsic  derivative  of  some  vector  bundle 
morphism  associated  with  the  DAE.  The  definition  we  now  give  is  an  equivalent  analytic 
translation  of  the  abstract  one  in  [RRh2],  specialized  to  the  DAE  (2.16)  -  (2.17).  The 
verification  of  this  equivalence  is  a  fairly  technical  exercise.  However,  the  method  used  in 
the  proof  of  [RIlh2,  Theorem  6.1],  dealing  with  a  special  case,  should  give  a  reliable  idea 
of  the  procedure  to  follow. 

Let  {t,  x)  he  given,  and  suppose  that 

dimker(5(t)Dj,G(t,x)|^^^^(,j  =  1, 

so  that  (t,  x)  G  by  condition  (2.12).  Evidently  we  have  rank  aco  ~ 

r  —  1,  whence 

dim[rge  Q(^)I>i(?(^ ^)|k„A(o]^  fl  Z  =  1. 

Now.  let  u  be  a  nonzero  vector  in  [rge  Q(t)DiG(t,  H  Z.  Eqmvalently,  u  G  Z  and 

[Q(t)D^G{t,x)]'^u  G  [ker  A{t)]^  =  rge  A(t)^,  so  that  there  is  a  unique  element  G  IR” 
such  that 

(2.19)  rje  rge  A{t),  A{t)'^ fj  =  [Q(t)D^G{t,x)]'^u. 

Remark  2.1:  Observe  that  the  condition  u  G  [rge  <3(t)DxG(i,  ,]■'■,  or  equivalently, 

u  G  ker(Q(t)Dj;G(f,x)[^^^^(,j)^,  in  no  way  implies  that  u  G  [rge  Q(t)DjG(t,x)]-‘-;  i.e., 
u  G  ker(Q(t)DxG(t,x))^,  as  could  inadvertantly  be  inferred.  Thus,  rj  in  (2.19)  may,  but 
need  not,  be  0.  □ 

These  preliminciries  lead  to  the  following  concept: 

Definition  2.3.  The  point  {t,x)  G  W  is  an  impasse  point  of  the  DAE  (2.1)  if 

(2.20)  dimker(3(t)DxG(t,x)|,„^(.,  =1, 

and  if 

(2.21)  {Q{t)DlGit,  x)(u)^  ft)  ^  0 
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holds  for  any  pair  of  nonzero  vectors 


(u,u)  e  x  {[rge  n  Zj. 

Remark  2,2:  (i)  Conditioa  (2.21)  is  unaiFected  by  the  choice  of  the  pedr  («,  ii)  since 
different  choices  amount  to  replacing  u  and  u  by  nonzero  scalar  multiples,  (ii)  Conditions 
(2.20)  and  (2.21)  are  also  independent  of  the  choice  of  the  space  Z  in  (2.3),  as  they  are 
equivalent  to  the  intrinsic  conditions  given  in  Definition  5.1  of  [RRh2].  □ 

The  definition  of  impasse  points  is  incomplete  without  the  notion  of  accessibility  based 
on  the  following  result: 

Proposition  2.2.  Let  (t,x)  be  an  impasse  point  of  (2.1)  (hence  (t,x)  G  W)  and  let  u  he 
as  in  Dehnition  2.3  and  dehned  by  (2.19).  Then 

(i)  the  vector 

(2.22)  i{D,(QG)(t,x),u},7j}  eRx  rge  A{t) 

is  nonzero  and  generates  the  orthogonal  complement  in  R  x  rge  A{t)  of  the  space  of  vectors 
of  the  form  (r,  A(f)/i)  with  (r,  h)  G  T(^t,x)W  C  R  x  R",  and 

(ii )  we  have 

(2.23)  (A(QG)(t,x),fi)  +  (G(i,x),rj}  ^  0. 

Proof,  (i)  By  contradiction,  suppose  that  t?  =  0,  i.e.  u  G  ker(Q(<)Dj;(j(t,  x))^  by  (2.19), 
and  suppose  that  (Dt(QG)(t,x),u}  =  0.  This  implies  that  u  G  [rge  D(QG)(t,x)]-^,  and 
hence  that  u  =  0  by  condition  (2.5)  since  u  £  Z.  This  contradicts  the  assumption  u  ^  0. 

It  is  easily  checked  that  the  vector  (2.22)  is  indeed  orthogonal  to  all  vectors  (r,  A(t)h) 
with  (t,  h)  G  and  it  follows  easily  from  condition  (2.21)  that  this  space  has 

dimension  r  —  1,  hence  codimension  1  in  R  X  rge  A(t).  This  proves  (i). 

(ii)  Once  again  by  contradiction,  suppose  that 

{Dt{QG){t,  x),  Z)  +  (G(t,  x),  n)  =  0, 
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so  that  the  vector  (l,G(t,x))  is  orthogonal  to  the  vector  (2.22).  As  (t,i)  G  W,  we  have 
Qit)G{t,x)  =  0,  i.e.  G(t,x)  €  rge  A{t),  and  it  follows  from  part  (i)  that  (l,G(t,x)) 
has  the  form  (r,A(t)p)  for  some  pair  (t,p)  €  Obviously,  r  =  1  and  hence 

(t,x,  l,p)  e  TW.  On  the  other  hand,  since  G{t,x)  =  A{t)p,  we  also  have  (f,x,  l,p)  G  M 
(see  (2.9)).  Thus,  (t,x,l,p)  G  TW  Cl  M,  whence  (f,x)  G  Tl{TW  D  M)  =  But 

then  Qit)D^G(t,x)i^^^^^^^  G  Gi:(ker  A(t),  Z)  by  condition  (2.12),  in  contradiction  with 
condition  (2.20).  □ 

From  condition  (2.21)  and  Proposition  2.2,  the  number 

(2.24)  {{D,iQG){t,xlu)  +  {G{t,x),p)){Qit)DlG(t,x)iu)\u) 

is  nonzero  when  (t,x)  is  an  impasse  point  of  the  DAE  (2.1).  We  shall  say  that  (t,  x)  is 
accessible  (resp.  inaccessible)  if  the  number  (2.24)  is  positive  (resp.  negative).  This  makes 
sense  because  of  the  following  result: 

Proposition  2.3.  If  (t,x)  G  W  is  an  impasse  point  of  (2.1),  the  sign  of  (2.24)  is  inde¬ 
pendent  of  the  complement  Z  of  rge  A(t),  and  independent  of  the  choices  of  u  and  u  in 
kerQ(t)D,G(f,x),,^^^^.,  \  {0}  and  [rge  ^  ^  ^  respectively.^ 

Proof.  .4s  various  arguments  are  involved  in  this  proof,  we  proceed  in  several  steps. 

(i)  For  fixed  Z,  the  expression  (2.24)  is  homogeneous  of  degree  2  in  w,  and  also  in  u 
since  fj  depends  linearly  upon  u,  so  that  its  sign  is  independent  of  the  choices  of  u  and  u. 
Moreover,  the  space  her  Q{t)DxG{t,x)\^^^ G  kerA(t)  :  DxG{t,x)h  G  rge  A(t)}  is 
independent  of  Z,  so  that  u  may  be  fixed  once  and  for  all. 

(ii)  The  expression  (2.24)  may  be  rewritten  in  a  way  making  no  use  of  the  derivative  DtQ{t): 
Since  (t,x)  G  W,  there  is  p  G  such  that  A{t)p  =  G(t,x),  whence  DtQit)G{t,x)  = 
—  Q{t)DtA{t)p  as  was  seen  in  the  proof  of  Proposition  2.1.  As  a  result,  we  have 

{DtiQG){t,x),u)  -  {Q{i)iDtGit,x)  -  DtAii)p),u). 

^Recall  that  rj  in  (2.24)  is  uniquely  determined  by  u  via  (2.19). 
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(iii)  Let  Zo  :=  [rge  and  call  (5o(*)  the  corresponding  (orthogonal  in  this  case) 

projection  operator.  It  is  straightforward  to  check  that  Q{t)'^  €  £(R")  is  also  a  projection 
operator  onto  Zq  (but  not  orthogonal  unless  Z  =  Zo).  As  a  result,  Q{t)'^u  e  Zq.  We 
claim  that  Q(t)^u  €  rge  the  choice  of  u,  we 

find  {u,Q(t)DxG{t,x)h)  =  0,  for  all  h  6  kerA(t),  i.e.  {Q{t)'^u,DxG{t,x)h)  =  0,  for  all 
h  e  ker  A(<).  As  Qo{t)  is  the  orthogonal  projection  onto  Zo,  Qo(t)  =  Qo{i)^  holds.  Hence, 
using  Q{t)'^u  =  Qoii)Qitfu,  we  get  {Q{tfu,Qo{t)D^G{t,x)h)  =  0,  for  aU  /i  €  ker  A(<), 
which  proves  the  claim.  At  this  stage,  we  have  that  uq  :=  Q{t)'^u  satisfies  the  condition 
5o  €  [rge  n  Zq.  Also,  from  (2.19),  rj  coincides  with  the  value  rjo 

obtained  by  choosing  Z  —  Zq  and  u  =  ug  in  the  first  place.  Using  (ii)  above  twice  (once 
with  Z,  once  with  Zq)  and  the  fact  that  Qo{t)  is  an  orthogonal  projector,  it  follows  that 
(2.24)  is  unchanged  when  replacing  u  by  Uq-  (hi  particular,  uq  ^  0  since  (2.24)  is  nonzero 
as  noticed  earlier.)  This  shows  that  in  (2.24)  we  may  replace  Z  by  Zo  upon  replacing  u 
by  uq.  But  from  (i),  «o  can  be  replaced  by  any  scalar  multiple  without  changing  the  sign 
of  (2.24).  Thus,  the  sign  of  (2.24)  is  independent  of  Z.  This  completes  the  proof.  □ 

Because  of  (2.20)  we  have  (t,r)  €  W  \  for  any  impasse  point  {t,x),  and,  hence, 
by  Theorem  2.1(i),  no  C’  solution  of  (2.1)  may  pass  through  {t,x).  For  this  reason,  we 
generalize  the  notion  of  a  solution  for  the  DAE  (2.1)  near  impcisse  points. 

Definition  2.4.  Let  the  DAE  (2.1 )  be  geometrically  nonsingular  of  index  1,  and  let  (t*,  z*) 
be  an  impasse  point  of  (2.1).  A  solution  of  (2.1)  satisfying  the  condition  x(t^)  =  z,  is  a 
continuous  function  z  :  J,  — >  R"  where  J*  =  +  T)  or  J*  =  (t.  —  T,t*]  for  some 

T  >  0,  such  that  z(i,)  =  z»  and  x  is  a  solution  of  (2.1)  in  J°  ~  J.^\  {t*}. 

Thus  solutions  of  (2.1)  in  the  sense  of  Definition  2.1  are  “one-sided”  and  need  not  satisfy 
A{t)x{t)  =  G(<,z(<))  for  i  =  i*.  For  the  proof  of  the  corresponding  existence  result,  given 
below,  we  refer  to  [RRh2]). 

Theorem  2.2.  Let  (<*,z,)  be  an  accessible  (resp.  inaccessible)  impasse  point  of  the 
geometrically  nonsingular  DAE  (2.1)  of  index  1.  There  are  exactly  two  solutions  x(t)  of 
(2.1)  in  the  sense  of  Dehnition  2.4  satisfying  the  conditions  x(t*)  =  z*,  and  both  are  dehned 
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in  J.  =  (t»  -  T,t^]  (resp.  [f„t,  +  T))  for  some  T  >  0.  This  result  remains  unaffected  by 
shrinking  T  >  0  and  both  solutions  are  actually  of  class  (7°°  in  J*  =  J*  \  {<»}.  Moreover, 
Hm  |ji:(t)||  =  oo.’^ 

This  theorem  justifies  the  terminology  “impasse  point”,  at  least  in  the  accessible  case 
since  the  solutions  cannot  be  continuously  extended  beyond  t, .  As  accessible  points  become 
inaccessible  and  vice-versa  upon  changing  time  evolution  (i.e.  changing  t  into  — <),  the 
terminology  is  justified  in  the  inaccessible  case  as  well.  Note  also  that  from  Theorems  2.1 
and  2.2.  impasse  points  lie  in  the  closure  of  in  W  despite  the  fact  that  such  a  property 
is  not  explicitly  incorporated  into  Definition  2.3.  Finally,  the  result  ^lim]jx(<)|l  =  oo  in 
Theorem  2.2  fully  justifies  dropping  the  requirement  that  A(t)x(f)  =  G{t,x{t))  for  t  —  t, 
in  Definition  2.4. 

Remark  2.3:  In  the  proof  of  Theorem  3.2  in  the  next  section,  we  shall  make  crucial  use 
of  a  property  more  precise  than  ^lim  ||i(f)||  =  oo  in  Theorem  2.2.  A  careful  examination 
of  the  proof  of  [R,  Theorem  5.1],  from  which  Theorem  2.2  is  derived  in  [RRh2],  reveals 
that  if  is  (say)  an  accessible  impasse  point  and  x  ;  (t*  —  T,  U]  —*  M"  denotes  either 

of  the  two  solutions  of  (2.1)  in  the  sense  of  Definition  2.4  and  satisfying  x(t*)  =  x„  then 
|[i(i)j|  =  0((t.  -f)"^/^)  for  i  near  This  shows  that  while  i(t)  blows  up  as  t  approaches 
t*,  nevertheless  x  €  —  T’,  f*))”.  □ 

Remark  2.4  (Autonomous  problems):  When  the  time  variable  t  does  not  enter  explicitly 
in  the  DAE  (2.1),  it  becomes  involved  rather  artificially  in  the  various  concepts  discussed 
earlier.  For  instance,  with  A(<)  =  .4  and  G{t,x)  =  G{x),  the  set  W  in  (2.6)  becomes 
Ty  =  Rx{x6lR”:  G{x)  €  rge  A}.  Here,  the  factor  R  is  needed  to  take  the  variable  t 
into  account,  but  the  only  “important”  factor  is  the  set  {x  G  R”  :  G(x)  €  rge  A}.  This 
suggests  changing  the  notation  for  autonomous  problems,  thus  eliminating  the  factor  R 
and  setting 

(2.25)  W  =  {x  G  M”  :  G{x)  G  rge  A}. 

^This  statement  differs  from  that  in  [RRh2,  Theorem  5.1]  where  typographical  errors  resulted  in  an 
exchange  of  the  roles  played  by  accessible  and  inaccessible  points. 
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Consistent  with  this  new  notation,  the  manifold  M  in  (2.9)  becomes 
M  =  {(a:,p)  e  R"  X  R”  :  Ap  -  G{x)  =  0}. 

In  this  case,  the  definition  of  the  set  of  consistent  points  given  in  (2.11)  remains 
adequate  (but,  of  course,  looses  its  artificial  time  component).  Likewise,  time  may  be 
dropped  from  the  definition  of  impasse  points  (Definition  2.3)  since  now  the  projection 
Q{t)  need  not  depend  upon  i.  Hence,  with  W  as  in  (2.25)  above,  we  may  and  shcill  refer  to 
r(€  W)  being  an  impasse  point  of  the  DAE  Ax  =  G{x)  (so  that  x  is  an  impasse  point  in  this 
“new”  sense  if  and  only  if  (t,  x)  is  an  impasse  point  for  every  t  €  R  in  the  “old”  sense).  The 
accessibility /inaccessibility  criterion  remains  based  upon  the  sign  of  the  quantity  (2.24), 
now  independent  of  t  and  hence  reducing  to  (G(r),^){QD^G(x)(u)^,u).  Naturally,  all  the 
results  discussed  earlier  have  obvious  einalogs  expressed  in  this  new  terminology.  □ 

3.  Discontinuous  Solutions  of  Semilinear  DAE’s. 

In  the  theory  of  explicit  ODE’s  x  =  /(t,  x)  with  smooth  enough  /,  it  is  rather  immaterial 
whether  x  should  be  viewed  as  the  classical  or  the  distribution  derivative  of  x,  except 
that  the  latter  point  of  view  introduces  a  few  extra  technicalities  since  f{t,x)  must  be 
unambiguously  defined  for  x  in  the  chosen  class  of  distributions.  In  particular,  viewing 
X  as  the  distribution  derivative  of  x  does  not  allow  for  the  existence  of  new  solutions  in 
the  class,  say,  of  piecewise  functions.  Indeed,  if  x  has  “jumps”,  these  get  multiplied  by 
Dirac  delta  distributions  in  x,  and  hence  must  vanish  for  i  to  equal  the  function  /(<,  x). 
It  follows  at  once  from  this  remark  that  piecewise  solutions  of  x  =  f{t,x)  are  just  its 
classical,  solutions.  We  shall  see  here  that  things  go  quite  differently  for  semilinear 
DAE’s. 

In  the  remainder  of  this  section,  we  shall  assume  once  and  for  all  that  the  DAE  (2.1)  is 
geometrically  nonsingular  of  index  1.  The  C*  solutions  of  (2.1)  about  consistent  points  as 
well  as  the  “one-sided”  solutions  of  (2.1)  about  impasse  points  in  the  sense  of  Definition 
2.4  will  be  referred  to  as  “classical”  solutions  of  (2.1).  The  solutions  we  shall  be  interested 
in  here  are  only  “piecewise  classical”  and  may  exhibit  jumps  at  one  or  several  points  of 
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their  domain  of  definition.  Naturally,  the  consideration  of  such  solutions  dictates  viewing 
i  as  a  generalized  derivative  of  x.  Not  surprisingly,  we  shall  choose  x  to  represent  the 
derivative  of  x  in  the  sense  of  distributions. 

If  J  =  (a,  6)  and  x  :  J  ^  R"  is  a  function  of  class  in  each  subinterval  (a,io]  and 
[to,  6)  for  some  to  €  J,  then  x  e  (I>'(  J))”  (distributions  in  J  with  values  in  R”)  and,  as  is 
well-known 

I  _  dx 

(3.1)  i  =  )^to +  -^, 

where  (resp.  xj )  =  lim  x(t)  (resp.  lim  x(f)),  is  the  Dirac  delta  distribution 

t— ►t(, 

at  to,  and  dx/dt  denotes  the  function  equal  to  the  usual  derivative  of  x  in  the  union 
(a,  to)  U  (to,i)  =  J  \  {to}- 

Remark  3.1:  Formula  (3.1)  need  not  be  true  if  x  is  in  (o,to)  and  (fo,&)  and  in 
(a,tol  and  [to,  6),  but  it  remains  valid  if  x  is  absolutely  continuous  in  (a,  to]  and  in  [to,&). 
This  generalization  will  be  crucial  to  the  proof  of  Theorem  3.2  later.  □ 

With  X  :  /  ^  R"  being  as  before  Remark  3.1,  the  frmction  G(t,  i)  also  defines  an 
element  of  (X>'(  J))”  in  the  obvious  way.  As  a  result,  A(t)x  -  G(t,  x)  is  a  distribution  on  J 
with  values  in  R”,  and  it  makes  sense  to  ask  whether  A{t)x  —  G(t,  x)  =  0  in  (P'(  J))”,  i.e. 
whether  x  solves  the  DAE  (2.1)  in  the  sense  of  distributions.  A  first  answer  is  given  next. 


Theorem  3.1.  Let  J  =  (a,  h)  and  let  (to,  Xq  )>  (^o,  )  G  be  given.  After  shrinking  the 

interval  J  about  to  if  necessary,  denote  by  x~  fresp.  x"^)  the  unique  solution  of  (2.1) 
in  J  satisfying  the  condition  x”(to)  =  x^  {resp.  x'^(io)  =  )>  "^bose  existence  follows 

from  Theorem  2.1,  and  set 

{x~(t)  ,  a  <  t  <  to, 
x+(t)  ,  to  <  t  <  &. 

Then,  x  €  0^(0,  to];  R”)  H  C’([to,  6);  R”)  and  x  solves  (2.1)  in  the  sense  of  distributions  if 
and  only  if 


(3.3) 


^0  “^0  €kerA(to). 
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Proof.  Trivial  from  formula  (3.1)  and  the  relation  A{t)^  —  G{t,x)  =  0  in  J  \  {to}-  Evi¬ 
dently,  condition  (3.3)  expresses  the  vanishing  of  the  coefficient  i4(to)(a^J  ~  of  ^<o 
the  expression  A(t)i  -  G(t,a:).  □ 

Interestingly,  if  (say)  (t^a;^)  G  W'^  is  fixed  in  Theorem  3.1,  condition  (3.3)  along  with 
the  requirement  (to,a^^)  €  imply  that,  in  general,  the  possible  values  for  form  a 
discrete  set.  Indeed,  since  diraker  j4(fo)  =  n  —  r,  condition  (3.3)  amounts  to  solving  a 
system  of  r  scalar  equations,  while  the  condition  (to,a:^)  G  W  (since  C  W)  holds  if 
and  only  if  <5(<o)G(to»a;J)  =  0,  i.e.  xj  solves  another  set  of  n  —  r  scalar  equations.  Thus, 
in  all,  x^  must  solve  a  set  of  n  scalar  equations,  usually  independent  of  one  another  and 
hence  having  only  isolated  solutions.  Note  that  by  openness  of  in  W  (see  Proposition 
2.1)  the  more  stringent  condition  (to,^^)  €  may  rule  out  some  of  the  solutions  x^  but 
does  not  place  any  further  limitation  upon  the  dimensionality  of  the  set  of  solutions,  so 
that  “generically”  this  set  should  remain  discrete.  Furthermore,  while  the  choice  i  J  =  Xq 
is  always  available,  the  nonlinear  nature  of  the  problem  makes  it  possible  for  solutions  with 
xj  7^  to  exist. 

Remark  3.2:  If  the  DAE  (2.1)  is  linear,  i.e.  of  the  form 

(3.4)  A{t)i  +  B{i)x  =  h{i), 

with  A,  B  and  b  of  class  C®®  so  as  to  fit  into  the  setting  of  this  paper,  then  jumps  do  not 
occur  and  the  only  solutions  of  (3.4)  in  the  sense  of  distributions  axe  the  classical  ones  (see 
[RRh4]).  Assuming  the  coefficients  A  and  B  axe  smooth,  discontinuous  solutions  of  (3.4) 
can  be  obtained  only  if  the  right-hand  side  b{t)  itself  is  discontinuous,  at  least  in  the  index 
1  case  considered  here.  Some  generalizations  will  be  discussed  in  Section  4.  In  contrast, 
nonlinearity  alone  is  responsible  for  the  existence  of  the  discontinuous  solutions  obtained 
in  Theorem  3.1.  This  is  well  illustrated  by  Example  3.1  below.  □ 

Example  3.1:  For  n  =  2  consider  the  autonomous  DAE  of  the  form  (2.1)  with 

(3.5)  A(0  =  A=(^  J,  G{t,x)  =  G{x)=l  ^  ),  x  G 

\0  0/  \Xi-fX2-x|/ 
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It  is  straightforward  to  check  that  (3.5)  is  geometrically  nonsingular  of  index  1  with,  in 
the  simplified  notation  for  autonomous  problems  discussed  in  Remark  2.4, 

The  points  =  (2\/3/9,  — \/3/3)  and  ^2  =  (— 2\/3/9,  \/3/3)  are  accessible  and  inaccessible 
impasse  points  for  (3.5),  respectively. 

Figure  3.1  gives  a  plot  of  W,  the  impeisse  points  ^2i  of  some  point  x*"  = 
(x)",^^'')  6  that  we  will  choose  as  initial  value  r”(0)  for  the  solution  x~  of  (3.5). 
This  solution  is  defined  imtil  it  reaches  the  accessible  point  ^1;  that  is,  its  interval  of  defini¬ 
tion  is  (—00, 2\/3/9  — ij").  The  arrows  in  Figure  3.1  represent  the  direction  of  evolution  of 
the  solutions  of  (3.5)  as  time  increases  (consistent  with  the  accessible/inaccessible  nature 
of  Cl  and  C2,  respectively). 


FIGURE  3.1 


.4ccording  to  Theorem  3.1,  jumps  can  only  occur  in  the  direction  of  the  null-space  of 
A,  i.e.,  in  the  “vertical”  direction.  As  a  result,  the  solution  x~  cannot  be  pieced  together 
with  another  solution  x"*"  before  Xj  (t)  reaches  the  interval  (— 2\/3/9, 2\/3/9);  that  is, 
before  t  6  X  :=  (-2v/3/9  -  x{”2y/Z/9  -  x^").  Indeed  for  to  <  -2i/3/9  ~  xj”  the  vertical 
line  xi  =  a;7(to)  intersects  W  only  at  the  point  x"(to);  but  as  soon  as  to  €  X,  this  line 
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intersects  W  at  x'"(to)  :=  and  two  other  points.  Either  one  can  be  chosen  as  initial 
condition  Xq  at  time  to  for  a  solution  of  (3.5),  and  then  formula  (3.2)  defines  a  solution 
of  (3.5)  in  the  sense  of  distributions  exhibiting  the  jump  —  x^  at  time  to-  Since  to  €  X 
was  arbitrary,  infinitely  many  (even  uncountably  many)  distinct  discontinuous  solutions 
may  be  obtained  by  this  process.  Furthermore,  after  a  first  jump  at  to,  the  solution  may 
jump  again  at  any  time  ti  6  X  with  tj  >  to  and  evidently  may  even  do  so  infinitely  many 
times,  in  infinitely  many  ways.  This  is  illustrated  in  Figure  3.2. 

It  is  important  to  observe  that  all  these  discontinuous  solutions  emanate  from  the  same 
point  x‘"  at  t  =  0.  In  other  words,  uniqueness  of  solutions  of  initial  value  problems 
associated  with  (3.2)  breaks  down  completely  in  the  class  of  discontinuous  solutions.  This 
rather  disturbing  remark  ~  obviously  not  limited  to  the  specific  example  (3.5)  -  motivated 
the  material  contained  in  Part  II  of  this  2crticle. 


A  pleasant  feature  about  discontinuous  solutions  is  that  they  provide  the  only  reasonable 
way  of  getting  out  of  the  dead  end  created  by  accessible  impasse  points.  Indeed,  since 
classical  solutions  cannot  be  continued  beyond  accessible  impasse  points  (Theorem  2.2) 
they  must  jump  to  proceed  further.  That  this  may  be  possible  for  solutions  of  (2.1)  in  the 
sense  of  distributions  is  shown  in  the  subsequent  generalization  of  Theorem  3.1. 
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Theorem  3.2.  Theorem  3.1  remains  valid  with  the  following  modifications  of  its  hy¬ 
potheses:  (to,  Xq)  is  an  accessible  impasse  point  of  the  DAE  (2.1),  x  is  either  of  the  two 
solutions  of  (2.1)  satisfying  x~(to)  =  Xq  (Theorem  2.2)  and  the  interval  J  —  (a,  6)  is  such 
that  x~  is  defined  in  (a,<o]- 

Proof.  The  arguments  of  the  proof  of  Theorem  3.1  can  be  repeated  verbatim,  but  the  use 
of  formula  (3.1)  must  be  justified,  which  is  done  via  Remarks  3.1  and  2.3.  □ 

The  proof  of  Theorem  3.2  given  above  is  deceptively  short,  as  all  the  technicalities  are 
contained  in  the  results  quoted  in  Remark  2.3.  There  is,  of  course,  also  an  analog  of 
Theorem  3.2  when  is  an  inaccessible  impasse  point,  and  one  when  (*0,2:^)  is  an 

accessible  impasse  point  and  an  inaccessible  impasse  point. 

Remark  3.3:  By  Theorem  3.2,  jumps  at  the  impasse  point  must  occur  in  the 

null-space  kerA(to).  This  result  agrees  with  Takens’  assumption  in  his  work  [T]  on  dis¬ 
continuous  solutions  for  special  cases  of  (2.1)  (autonomous  with  a  gradient  structure)  and 
also  with  the  work  of  Sastry  and  Desoer  [SDe].  But  neither  Takens  nor  Sastry  and  Desoer 
ever  consider  distribution  solutions,  and  they  use  other  arguments  to  justify  their  concepts 
of  solution.  Also,  in  these  approaches,  jumps  are  allowed  to  occur  only  at  singularities. 
(In  [SDej,  it  is  observed  that  jumps  could  also  occur  at  other  points,  but  this  possibility  is 
next  ruled  out  by  the  introduction  of  extra  assumptions.)  □ 

It  should  be  pointed  out  that  Theorem  3.2  does  not  state  that  jumps  always  allow  the 
solutions  to  be  continued  beyond  impasse  points  as  is  obviously  the  case  for  Example  3.1. 
For  example,  Figure  3.3  below  relates  to  the  DAE 

fxi  =  l, 

\o  =  xi  +X2, 

for  which  A  is  as  in  (3.5)  and  W"  =  {(-r^,  2:2);  X2  G  R}  (using  once  again  the  “autonomous” 
notation  of  Remark  2.4).  In  this  case,  jumps  from  impasse  points  cannot  occur.  In  fact, 
x~  =  (0,0),  is  the  only  (accessible)  impasse  point  and  there  is  no  point  G  such  that 
2;+  _  2r^  €  ker  A,  and  hence  no  point  to  which  a  solution  reaching  x^  can  jump. 
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FIGURE  3.3 


4.  Forced  Discontinuities  and  Inconsistent  Initial  Conditions. 

The  discontinuous  solutions  discussed  in  the  previous  section  axe  in  some  sense  “self- 
generated”,  as  they  exist  from  only  the  combined  effect  of  the  nonlinearity  of  G{t,x)  with 
respect  to  the  x  variable  and  the  weaker  requirement  that  (2.1)  be  imderstood  in  the  sense 
of  distributions.  As  noted  in  Remark  3.2,  such  self-generated  discontinuities  do  not  exist 
in  linear  problems. 

A  different  situation,  also  frequently  encountered  in  practical  applications,  is  that  the 
DAE  governing  the  system  suddenly  changes  at  a  given  time  to-  Actually,  this  situation 
always  presents  itself  when  the  DAE  (2.1)  begins  to  govern  the  system  at  time  to,  and 
the  state  variable  has  evolved  in  an  unrelated  way  for  t  <  to.  Mathematically,  we  may 
assume  that  this  “past”  history  is  accounted  for  by  a  known  function  x~{i)  for  t  <  to  and 
that  lim  x~(t)  :=  Xq-  Evidently,  is  a  natural  initial  condition  to  associate  with  the 
DAE  (2.1)  if  the  latter  is  to  describe  the  state  of  the  system  for  t  >  to-  The  only  problem 
of  course  is  that  there  is  no  reason  why  (to»a:^)  should  satisfy  the  consistency  condition 
(to  j  a:o  )  ^  which  is  necessary  for  (2.1)  to  have  a  classical  solution  satisfying  i(to)  =  x^. 
This  directly  leads  to  the  well  known  problem  of  inconsistent  initial  conditions.  This 


problem  is  unsolvable  in  the  framework  of  classical  solutions,  and  has  a  straightforward 
answer  when  discontinuities  are  permitted  and  solutions  are  understood  in  the  sense  of 
distributions. 

Indeed,  the  problem  formulated  above  is  one  of  extending  a  known  function  x~{t)  into 
a  solution  of  (2.1).  Let  b{t)  be  the  function 


(4.1) 


6(t)  = 


0  for  to  <  ^  <  b, 


for  G  <  t  <  to 


where  dx  /dt  denotes  the  classical  derivative  of  x  (assuming  x  of  class  in  (a,  to])  and 
consider  the  problem  of  finding  a  function  x  of  class  in  (a,  to]  and  in  [fp,  b),  satisfying 


(4.2) 


{A(t)d:  -  G(t,  x)  =  b(t)  in  (V'(J))’^, 


It  is  easily  checked  that  x  solves  (4.2)  if  and  only  if  x(t)  =  x  (t)  for  t  <  to,  and  a;''' 

solves 


U(t)-^-G(t,x+)  =  0 
\  x+(to)  =  x+. 


with  x^  —  x^  €  ker.4(to)  (so  that  no  Dirac  delta  distribution  appears  in  (4.2)).  As  in 
the  previous  section,  we  thus  have  that  x^  is  determined  from  x^  by  the  two  conditions 
Xo  -Xq  e  ker.4(to)  and  Q(to)(?(to,  x^)  =  0,  the  latter  for  consistency  of  (to,3r^)  with  the 
DAE  (2.1)  (i.e.  (4.3)).  By  nonlinearity  of  G,  x^  is  usually  non-unique.  On  the  other  hand, 
from  the  fact  that  the  conditions  x^  —  x^  €  ker  A(to)  and  (9(to)G(fo,3:^)  =  0  represent 
a  system  of  n  equations  in  n  scalar  unknowns  for  x^,  only  isolated  solutions  should  be 
expected.  These  considerations  provide  a  simple  satisfactory  answer  to  the  problem  of 
inconsistent  initial  conditions  (although  not  quite  complete  since  xf  need  not  be  imiquely 
determined).  Further  limitations  on  the  choice  of  xj  will  be  dictated  by  the  results  in  Part 
II  of  this  article. 


22 


References 


[CG]  Campbell,  S.L.  and  Griepentrog,  E.,  Solvability  of  General  Differential-Algebraic  Eqnaiions, 
SIAM  J.  Sci.  Comp,  (to  appear). 

[CD]  Chua,  L.  O.  and  Deng,  A.-C.,  Impasse  Points.  Part  I:  Numerical  Aspects,  Int.  J.  Circ.  Th.  and 
Appl.  17  (1989),  213-235. 

[R]  Rabier,  P.J.,  Implicit  Differential  Equations  Near  a  Singular  Point,  J.  Math.  Anal.  Appl.  44 
(1989),  425-449. 

[RRhl]  Rabier,  RJ.  and  Rheinboldt,  W.C.,  A  Geometric  TVeatment  of  Implicit  Differential- Algebraic 
Equations,  3.  Diff.  Equations  109  (1994),  110-146. 

{RRh2]  Rabier,  P.J.  and  Rheinboldt,  W.C.,  On  Impasse  Points  of  Quasilinear  Differential-Algebraic 
Equations,  J.  Math.  Anal.  Appl.  181  (1994),  429-454. 

[RRh3]  Rabier,  P.J.  and  Rheinboldt,  W.C.,  On  the  Gomputation  of  Impasse  Points  of  Quasilinear 
Differential-Algebraic  Equations,  Math.  Comp.  62  (1994),  133-154. 

[RRh4]  Rabier,  P.J.  and  Rheinboldt,  W.C.,  Classical  and  Generalised  Solutions  of  Time- Dependent  Lin¬ 
ear  DAE’s,  Inst,  for  Comp.  Math,  and  Appl.,  Univ.  of  Pittsburgh.  Tech.  Rept.  TR-ICMA-183, 
October  1993,  Linear  Algebra  and  Appl.,  (to  appear). 

[RRh5]  Rabier,  P.J.  and  Rheinboldt,  W.C.,  Time- Dependent  Linear  DAE’s  with  Discontinuous  Inputs, 
Inst,  for  Comp.  Math,  and  Appl.,  Univ.  of  Pittsburgh,  Tech.  Rept.  TR-ICMA-186,  March  1994, 
Linear  Algebra  and  Appl.,  (to  appear). 

[SDe]  Sastry,  S.  S.,  and  Desoer,  C.A.,  Jump  Behavior  of  Circuits  and  Systems,  IEEE,  Trans.  Circ.  and 
Syst.  28  (1981),  1109-1124. 

[T]  Takens,  F.,  Constrained  Equations:  A  Study  of  Implicit  Differential  Equations  and  Their  Discon¬ 
tinuous  Solutions,  Lecture  Notes  in  Math.  Vol  525,  Springer- Verlag,  New  York,  1976,  pp.  143-234. 


23 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  WO.  0704-0788 


J  OATES  COVERED 

REPORT 

4.  TITLE  AND  SUBTITLE 

DISCONTINUOUS  SOLUTIONS  OF  SEMI  LI  NEAR 
DIFFERENTIAL-ALGEBRAIC  EQUATIONS 

PART  1:  DISTRIBUTION  SOLUTIONS  .  - 

5.  FUNDING  NUMBERS 

ONR-N-00014-90-J-1025 

NSF-CCR-920348S 

6.  AUTHOR(S) 

Patrick  J.  Rabier 

Werner  C.  Rheinboldt 

7.  PERFORMING  ORGANIZATION  NAME(5)  AND  ADDRESS(ES) 

Department  of  Mathematics  and  Statistics 

University  of  Pittsburgh 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

ONR 

NSF 

10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

11.  SUPPLEMENTARY  NOTES 

^  12a.  DISTRIBUTION /AVAILABILITY  STATEMENT 

Approved  for  public  release:  distribution  unlimited 

12b.  DISTRIBUTION  CODE 

i  13  ABSTRACT  (Maximum  200  word!)  ■  i  u  • 

j  There  is  strong  physical  evidence  that  a  full  treatment  of  differential-algebraic 
\  equations  should  be  incorporate  solutions  with  jump  discontinuities.  It  is  shown 
i  here  that  for  semi  I  inear  problems  the  setting  of  distributions  allows  for  the 

development  of  a  theory  where  indeed  such  discontinuities  may  occur.  This  approach 
also  settles  the  problem  of  inconsistent  initial  conditions  in  a  very  simple  way. 

On  the  other  hand,  new  issues  arise  as  not  only  uniqueness,  but  even  countability 
of  the  number  of  solutions  of  initial  value  problems  may  now  be  lost.  A  physically 
motivated  but  purely  mathematical  selection  procedure  to  overcome  this  difficulty 
is  discussed  in  Part  II  of  this  paper. 


I 

! 

j 


3 

4 


i 

■; 

J 


I 


14.  SUBJECT  TERMS 

Differential-algebraic  equations,  semi  linear,  discontinuous 
solutions,  distribution  soiutions 

15.  NUMBER  OF  PAGES 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

OF  REPORT 

unc lassi f ied 

18.  SECURITY  CLASSIFICATION 

OF  THIS  PAGE 

unclass i f i ed 

19.  security  CLASSIFICATION 
OF  ABSTRACT 

unclassi f ied 

20.  LIMITATION  OF  ABSTRACT 

